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Abstract 

We propose a novel quasiparticle interpretation of the equation of state of deconhned 
QCD at finite temperature. Using appropriate thermal masses, we introduce a phenomeno- 
logical parametrization of the onset of confinement in the vicinity of the predicted phase 
transition. Lattice results of the energy density, the pressure and the interaction measure 
of pure SU (3) gauge theory are excellently reproduced. We find a relationship between 
the thermal energy density of the Yang-Mills vacuum and the chromomagnetic condensate 
(B^)t- Finally, an extension to QCD with dynamical quarks is discussed. Good agreement 
with lattice data for 2, 2+1 and 3 fiavour QCD is obtained. We also present the QCD 
equation of state for realistic quark masses. 
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1 Introduction 



There are convincing arguments that QCD exhibits a transition from a confined hadronic phase 
to a deconfined partonic phase, the quark-gluon plasma (QGP), at a temperature of Tq ~ 170 
MeV (for two hght quark flavours) [|l|, ^. A central quantity of matter in thermal equilibrium 
is the Helmholtz free energy from which the pressure p, energy density e and entropy density 
s can be derived. These entities are important for the description of, e.g., ultra-relativistic 
heavy-ion collisions at the CERN SPS and RHIC, or the evolution of the universe after the 
first 10~^ — 10~^ seconds. The challenge lies therefore in the derivation of the equation of 
state (EoS) of hot, deconfined QCD from first principles. Perturbative results are available 
up to order 0{g^) Q. However, for temperatures of interest in the experimentally accessible 
region (a few times Tc), the strong coupling constant is presumably large: gs — 1.5 — 2. The 
perturbative expansion in powers of gs shows bad convergence already for much smaller values 
of the coupling. Furthermore, in the vicinity of a phase transition, perturbative methods are in 
general not expected to be applicable. Non-perturbative methods such as lattice QCD calcula- 
tions become mandatory. From these numerical simulations the EoS of a pure gluon plasma is 
known to high accuracy |Q, |5| , and there are first estimates for the continuum EoS of systems 
including quarks, albeit still with unphysically large masses [^. 



Various interpretations of the lattice data have been attempted in terms of physical quantities, 
most prominently as the EoS of a gas of non-interacting, massive quark and gluon quasiparti- 
cles. Their thermally generated masses are based on perturbative calculations carried out in 
the hard thermal loop (HTL) scheme 0, More recently, the QGP has also been described 

in terms of a condensate of Wilson lines |1C]. Since all microscopic dynamics has been 
integrated out in the EoS, there exists no unique interpretation of the lattice data, and one 
must resort to additional information in order to further restrict the setup of such models. The 
phenomenological models of the QGP as an ideal gas of massive quarks and gluons have found 
support from resummed perturbation theory |11| for temperatures T > 3 Tc- However, it is 
astonishing that the lattice EoS is described in this model even close to Tq where one would 
expect non-perturbative dynamics to enter. Currently there exists not even a qualitative mi- 
croscopic explanation for this behaviour, and the model has difficulties explaining the dropping 
of the thermal gluon screening mass in the vicinity of the phase transition. 



In this work we extend the quasiparticle approach. Our main new ingredient, as compared 
to previous work, is a phenomenological parametrization of (de) confinement, supplemented by 
thermal quasiparticle masses compatible with lattice results. In section II, we shortly review the 
quasiparticle formalism. Our model is set up, and results for various observables are compared 
with lattice data of continuum extrapolated SU(3) gauge theory. In section HI we give argu- 
ments how an extension to a theory including quarks might work, show how the existing lattice 
data are fitted and estimate the QCD EoS for realistic quark masses. Section IV summarizes 
and gives an outlook. 



2 Gluonic quasiparticles 
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2.1 Basic quasiparticle model 



In this section, we consider an SU{Nc) gluon plasma {Nc = 3) at finite temperature and 
vanishing chemical potential. The use of a quasiparticle model in QCD is based on the obser- 
vation that in a strongly interacting system, the complex dynamics often rearranges itself in 
such a way that gross features of the physics can be described in terms of appropriate effective 
degrees of freedom. From asymptotic freedom, we expect that at very high temperatures the 
plasma consists of quasifree gluons. As long as the spectral function of the thermal excitations 
at lower temperatures resembles qualitatively this asymptotic form, a gluonic quasiparticle de- 
scription is expected to be applicable. The dispersion equation for transverse gluons reads 
a;^ — A;^ — Il^{LO,k) = 0. Here, k = \k\, and IIj is the transverse part of the thermal gluon 
self-energy. If, for thermal momenta uj,k T, the momentum-dependence of 11* is weak and 
its imaginary part small, gluon quasiparticles will propagate mainly on-shell with the dispersion 
relation 

u\k)c^k^+ml{T), (1) 

where mg{T) acts as an effective mass generated dynamically by the interaction of the gluons 
with the heat bath background. Since the existence of a preferred frame of reference breaks 
Lorentz invariance, new partonic excitations, longitudinal gluonic plasmons, are also present in 
the plasma. However, their spectral strengths are exponentially suppressed for hard momenta 
and large temperatures, so gluons are expected to retain their Vg = 2{Nq — 1) degrees of freedom 
despite their masses. 



For homogeneous systems of large volume V, the Helmholtz free energy F is related to the 
pressure p by F{T,V) = —p{T)V. In the present framework of a gas of quasiparticles, its 
explicit expression reads 

P{T) = dk fB{E,)— - B{T), (2) 

where is the gluon degeneracy factor, Ef^ = ^Jl^~-\-m^{T) and 

fsiEk) = , , ^ -■ (3) 

exp{Ek/T) - 1 

The energy density e and the entropy density s take the form 

e{T) = ^J^ dk k^fBiEk) Ek + B{T) (4) 

and 

The function B{T) is introduced to act as a background field. It is necessary in order to 
maintain thermodynamic consistency: eqs.(^), and (|5|) have to satisfy the Gibbs-Duhem 
relation e + p = sT, and s(T) is related to p(T) via 
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B{T) basically compensates the additional T-derivatives from the temperature-dependent masses 
in p and is thus not an independent quantity. Since B{T) adds to the energy density of the 
quasiparticles in eq.@, it can be interpreted as the thermal vacuum energy density. The 
entropy density, as a measure of phase space occupation, is unaffected by B(T). 



2.2 The HTL model 

In previous work ^, the HTL perturbative expression 



mg{T) = J ^ 9s{T)T with g'.iT) = (7) 

has been used to model the thermal gluon mass^ Here, mg{T) follows from the transverse part 
HJ" of the thermal gluon polarization tensor in the limit k S> QsT and it is gauge-independent 
to this order. Phenomenology enters in the effective coupling constant gs(T) with the two 
fit parameters Tg and A. By identifying the Landau pole of the effective coupling with a 
temperature close to Tc, the effective gluon mass (^) becomes very heavy in the vicinity of 
the phase transition, and s, p and e drop abruptly to almost zero. The resulting EoS of the 
HTL model is then in good agreement with lattice data over a wide temperature range between 
Tc and 5 Tc Q. For temperatures T > 3 Tc, HTL resummed perturbation theory indeed 
supports the picture of the QGP as a gas of weakly interacting, massive quasiparticles. However, 
there is no a priori reason to expect that the one-loop expressions for gs{T) and H^ can be 
extrapolated to values of the coupling as large as gs{Tc) ~ 2. Furthermore, close to a phase 
transition, the reliability of perturbative calculations is questionable. In the light of these facts, 
it is astonishing how well the HTL model works even in the vicinity of Tq- 



Numerical simulations suggest that the deconfinement transition in Yang- Mills theories is second 
order for Nq = 2 |16| and weakly first order for the physical case Nq = 3 From the 

general theory of critical phenomena, it is expected that the correlation length ^(T), which 
is proportional to the inverse of the gluonic screening Debye mass m^), grows when Tc is 
approached from above (we recall that rriD measures the exponential decay of the static gluon 
field correlator {A'Q{r)A\^{{)))T oc (5"'' exp(— m£)(T)|f|)/|r|). For three colours, the mass gap does 
not vanish at Tc, so i{Tc) remains large, but finite. This behaviour is indeed seen in lattice 
calculations ||l5| : m/j drops by a factor of ten when going down from 2 Tc to Tc (see figure 
0). 



In HTL perturbation theory, m£,{T) and mg{T) are connected by the simple relation: 

mo = V2 rrig. (8) 

A scenario with heavy masses rUg would then imply small correlation lengths close to Tc- It 
is therefore not clear how a decreasing gluonic Debye mass can be matched to the heavy, non- 
interacting quasiparticles of the HTL model. Of course, in a more general non-perturbative 
framework, rUg and m^i, although both arise from the same polarization tensor, will not be 
related a priori by a simple constant, as in eq.(^). 

''In the following, we refer to the quasiparticle model with the gluon mass defined in eq. (0) as the HTL (hard 
thermal loop) model. 
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2.3 Deconfinement and the Quasiparticle model 



We believe that the failure of simple quasiparticle models to exhibit the correct temperature de- 
pendence of the Debye mass can be traced back to the fact that the picture of a non-interacting 
gas is not appropriate close to Tq because the driving force of the transition, the confinement 
process, is not taken into account. Below Tc the relevant degrees of freedom in a pure SU{2>) 
gauge theory are heavy, colour singlet glueballs. Approaching Tc, deconfinement sets in and the 
gluons are liberated, followed by a sudden increase in entropy and energy density. Conversely, 
when approaching the phase transition from above, the decrease in the thermodynamic quan- 
tities is not caused by masses becoming heavier and heavier, instead the number of thermally 
active degrees of freedom is reduced due to the onset of confinement. As T comes closer to 
Tc, an increasing number of gluons gets trapped in glueballs which disappear from the thermal 
spectrum: since mcB ^ 1-5 GeV and Tc ~ 270 MeV (for pure gauge theory), glueballs are 
simply too heavy to become thermally excited in the temperature range under consideration 
(up to about 5 Tc). Of course, glueball masses may also change with temperature. However, 
since the contribution of glueballs to e, p and s is negligible below Tc (as evident from lattice 
data), this change is presumably small. The important fact in our opinion is the following: 
while the confinement mechanism as such is still not understood, it is not necessary to know 
it in detail since we consider a statistical system. All confinement does on a large scale is to 
cut down the number of thermally active gluons as the temperature is lowered. The question 
is whether this effect of confinement can be reconciled somehow with the quasiparticle picture. 
We will show in the following that it is indeed possible in a simple, phenomenological way if we 
allow for an effective, temperature-dependent number of degrees of freedom Vg[T). 

Let us assume that the thermal gluon mass mg{T) does not increase as Tc is approached, but 
instead follows roughly the behaviour of the Debye mass, i.e. it decreases. Its detailed T- 
dependence is not important for the discussion at the moment, but it will be examined in more 
detail in the next section. Consider now the entropy of a gas of massive gluons along eq.(|5|) 
with such a dropping effective gluon mass. The result for s{T) will clearly overshoot the lattice 
entropy because light masses near Tc lead to an increase in s. However, since the entropy is a 
measure for the number of active degrees of freedom, the difference may be accounted for by the 
aforementioned confinement process as it develops when the temperature is lowered toward Tc. 
This effect can be included in the quasiparticle picture by modifying the number of effective 
degrees of freedom by a temperature-dependent confinement factor C{T): 




(9) 




(10) 





and 



s{T) 




oo 



dk e [C{T)fB{Ek)] 



^k^ + mliT) 
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In essence, the factor C(T) represents a statistical parametrization of confinement. Its explicit 
functional form can be obtained simply as the ratio of the lattice entropy and the entropy ( |l2|) 
calculated with a dropping input gluon mass mg{T). Qualitatively, we expect C(T ^ Tc) ~ 1 at 
high temperatures where the deviation from the Stefan-Boltzmann limit of e, s and p, as seen on 
the lattice, is caused solely by the thermal masses rUg. As the phase transition is approached 
from above, the number of thermally active degrees of freedom decreases and consequently, 
C{T) becomes less than one. Finally, the entropy below Tq is small, but non-zero, and we 
can estimate C{Tc) ~ 0.2 from lattice data. We would also expect that C{T) is a smooth, 
monotonously increasing function with T, following the behaviour of the entropy density. 

While the lattice entropy is fitted by construction, the crucial test of the model lies now in 
the reproduction of the energy density e and the pressure p which follow unambiguously from 
eqs.(pl]) and (|l^). Owing to eq.(^, the background field B{T) depends on mg{T) and C(T) 
through: 

B{T) = Bi{T) + B2iT) + Bo, where (13) 



T 



A 



Tc 
T 



B2{T) = J dr C{t)^^ J dk fB{Ek) 

Tc 



Setting C(T) = 1 va. B2 and i?i = 0, one recovers the HTL model expression for B{T). The 
integration constant Bq is chosen such that the gluonic pressure equals the very small glueball 
pressure pgb (which is taken from the lattice) at Tq, according to Gibbs' condition p^^^on = Pgb- 

Note that the proposed model for the EoS has no free fit parameters in the sense that once 
the non-perturbative T-dependence of the thermal gluon mass is fixed, C(T) follows from the 
ratio of the entropy density of lattice QCD and the entropy density calculated with mg{T). 
B{T) is fixed up to an integration constant, which is obtained from Gibbs' condition. Then 
the energy density and the pressure are uniquely determined. Since the pressure is related to 
the partition function Z by p{T) = —F{T,V)/V = T/VlnZ, the full thermodynamics of the 
system is known. 

2.4 Thermal masses 



We must now specify our input thermal quasiparticle mass mg{T) in eqs.(10 - 13). As mentioned 
earlier, mg{T) is to be identified with the transverse part Il*{uj,k;T) of the gluon polarization 
tensor at u;, /c ~ T (see eq.(|l|)). Evaluating mg[T) requires a detailed non-perturbative analy- 
sis of the gluonic two-point correlation function which lattice calculations could, in principle, 
provide. In practice this information does not (yet) exist, so we have to rely on a model. 

Suppose we still keep the basic form of eq.(0), 

mg{T) = G{T)T, (14) 
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but assume that the dimensionless effective coupHng G{T) shows approximate critical behaviour 
at T close to Tq'- 

G{T)^Go(l-^y , (15) 



with some characteristic exponent (3 and a constant Go -The assumption ( |15| ) implies that the 
thermal mass behaves as m„ (T) ~ (T - Tc)^ close to Tc- Asymptotically at T > Tc, G{T) 
should match the HTL perturbative form as in (0). In practice we can choose this matching 
point, for instance, at Tm = 3 Tc- This fixes Go ^JNc/& gs{Tm) — 1-3. A similar value for 



Tfn below which an explicit HTL resummation is expected to fail, was obtained in ref.|12| 



The quantity for which lattice information does exist is the Debye screening mass mD{T) that 
is related to the longitudinal part 11^(0;, /c;r) of the polarization tensor. When defined as 
m|) = n^*(0, fc^ = —mj-)) ||l3|, the result turns out to be gauge-independent for a wide class of 
gauges to arbitrary order in perturbation theory (unlike the situation for lo = 0, k ^ 
where H^* is not gauge-invariant at next-to-leading order). While there is no a priori reason why 
rUg and m£) should still be related non-perturbatively as they are in perturbation theory (c/. 
eq.(P)), it is nevertheless instructive to recall what is known about the temperature dependence 
of the Debye screening mass above Tc- Explicit values for have been extracted from lattice 
calculations of (colour-averaged) heavy quark potentials V{R,T) |]15[ by the ansatz 



^ (16) 



Perturbation theory predicts d = 2 and n{T) = 2m ^{T)- As elucidated in |15], the potential 
may be better reproduced in terms of a mixture of one- and two-gluon exchange since the 
observed behaviour close to Tc favours values of d ~ 1.5 in eq.(|l6|). It is now interesting to 
observe that the lattice result for /i(T) can be parametrized very well by 

/x(r) const. -T ('[1 + 5]-^')'' (17) 



with /5 ~ 0.1 and a small gap at T = Tc introduced by 5 10 . The form of eq.(|l7[) is indeed 
reminiscent of approximate critical behaviour. 

Let us then assume that the proportionality (^) between the screening mass rriD and the thermal 
gluon mass rUg remains at least qualitatively valid in the vicinity of the phase transition, i.e. 
that the exponent /? in the characteristic (1 — Tc /T)^ behaviour of both rriD and nig is roughly 
the same. As it turns out, this is not a serious assumption: we have checked that, as long as 
mg(T) and mj:i{T) just have similar trends in their T-evolution close to Tc, our results are not 
sensitive to the detailed quantitative behaviour of the quasiparticle mass. 

Guided by these considerations, the thermal gluon mass mg is thus parametrized as 

mg{T) = GoT(^[l + 6]-^y . (18) 



where we allow for a small mass gap at T = Tc, as indicated by the lattice results for moiT) |15| . 
The small correction 6 <^ 1 encodes this deviation. Finally, Go is determined by the asymptotic 
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Go 


6 


(3 


Co 


Sc 




Set A 


1.35 


10"^ 


0.2 


1.24 


0.0029 


0.34 


Set B 


1.30 


10-6 


0.1 


1.25 


0.0026 


0.31 


Set C 


1.30 


10-7 


0.05 


1.27 


0.0021 


0.30 



Table 1: Parametrizations for the thermal gluon mass mg{T)/T and the corresponding confine- 
ment factor C(T). 



value of the thermal mass, chosen such that the lattice mass and the HTL perturbative result 
from 1^ coincide at T ~ 3 Tc, as mentioned before. In order to account for uncertainties and 
the approximate nature of relation (p!8|), we have investigated a range of values for Go, S and /? 
which can be found in table ||. The upper limit of the range is labelled Set A, an intermediate 
parameter set Set B and the lower limit Set C. Its plots and the lattice data points for m£,[T) 
are displayed in figure ||. 

A decreasing effective coupling strength G(T) as Tq is approached from above, seems at first 
sight counterintuitive: One would expect that, at a scale T ~ ^qcd-, 'infrared slavery' sets in, 
accompanied by an increasing QCD coupling gg. However, it should be borne in mind that this 
expectation is based on a perturbative result extrapolated to large couplings, neglecting non- 
perturbative effects. A heuristic argument to make the dropping effective coupling plausible 
goes as follows. Since we are in a strong coupling regime, the interactions between gluons 
cannot be described in terms of single gluon exchange, instead they are dominated by (non- 
perturbative) multi-gluon dynamics. As the temperature is lowered, more and more gluons 
become confined and form heavy glueballs, as outlined earlier. The effective glueball exchange 
interaction between gluons reduces approximately to a local four-point interaction proportional 
to I/tuq^. The total interaction can be interpreted as a superposition of multi-gluon and 
(weak) glueball exchange. Obviously, the more glueballs are formed, the weaker becomes this 
interaction. The coupling G(T) in eqs.(|l^) and (18) reflects an interaction between bare gluons 



from the heat bath on length scales l/mcB and larger that turns these bare gluons into massive, 
weakly interacting quasiparticles on length scales of order 1/T. 



2.5 Results for SU(3) gauge theory 

Now that the temperature behaviour of mg{T) is given, we can explicitly calculate the entropy 
density (|T^). Dividing the lattice entropy density by the result of this calculation, we obtain 
the T-dependence of the confinement factor G(T). A very good fit of the resulting curves again 
exhibits an approximate critical power-law behaviour: 

G(r) = Go([l + 5c]-^)^\ (19) 

which is a non-trivial result. The corresponding parameters of G(T) for the different mass 
parametrizations, sets A, B and C of (0), can also be found in table 0. Their plots are shown 
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Figure 1: The thermal gluon mass rrig/T. The grey band shows the parameter range of table ^ 
for eq.(18). Set A marks the upper limit and set C the lower, the intermediate set B is displayed 
as a solid line. The arrows indicate the mass gap mg{Tc) /Tc of the different sets at the critical 
temperature. Symbols display ijl{T)/2 of eq.([T^), with d = 1.5 fixed, for lattice configurations 
of different temporal extent N^- [l^. For comparison the HTL perturbative mass m^^^/T of 
eq.(0) is also plotted (dashed line), using the parameters of ref.[^. 

in figure |2[ C (T) is obviously only weakly sensitive to variations of the mass parameters within 
a broad band, so the further discussion will be based on Set B. In the following, the quasiparticle 
model with the gluon mass (18) and the confinement factor ( |l9| ) is referred to as 'confinement 
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model'. For the integration constant Bq appearing in eq.(p!3|) we find 0.30 ~ (200 MeV) 
This value is about a factor of 2 larger than in the HTL model and remarkably close to the value 
of the bag constant at T = 0, a welcome feature. As mentioned, we expect the confinement 
effect to be negligible for T > 3 Tq where the HTL quasiparticle model sets in. From there 
on, C(r) — > 1. The actual deviation of C{T) from 1 for large temperatures has two reasons: 
first, even in the HTL quasiparticle model a gluon degeneracy of = 16 does not describe the 
data, instead a value larger by about 10% is necessary to account for residual sub-leading effects 
not captured by the model. Second, the behaviour of the gluon mass for larger T is certainly 
oversimplified since the parametrization rrig oc (T — Tq)^ is expected to be valid only in the 
vicinity of Tq- Its value overestimates the HTL perturbative result for T > 3 Tc by some 5%, 
hence the thermodynamical potentials are slightly smaller than in the HTL model. One should 
instead apply a smooth interpolation between the HTL perturbative mass and the mg(T) we 
used for temperatures close to Tc- This would in turn yield a more complicated expression for 
C(T), but it is in principle straightforward. 

In figure ^ we compare results of the confinement model to continuum extrapolated SU (3) 
lattice data B. Obviously, the thermodynamic quantities are very well described even close to 
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Set B 



T/T, 



Figure 2: The confinement factor C(T) as a function of temperature. The grey band shows the 
range for the corresponding mass parametrizations of table |l|. The sohd fine is obtained from 
set B. 

Tc- The shght deviations in the region ~ 5 Tc arise from our simple parametrization of C{T). 
We want to stress again that the entropy density is, by construction, always fitted. The highly 
non-trivial test of the proposed model lies in the reproduction of e and p. 

A quantity that is sensitive to the finer details of the model is the trace of the energy- momentum 
tensor, T^^ = e — 3p, which is compared to data from a 32^ x 8 lattice in figure The interaction 
measure 

A(T) = (e - 3p)/T' (20) 
is connected, via the QCD trace anomaly, to the temperature dependent gluon condensate: 

r^A(r) = {g')t=o - {g')t. (21) 

{G')t=^-^{G%')t and {G')t=o = (2.5 ± 1.0)Tt- (22) 

Again, excellent agreement over the whole temperature range is observed. The confinement 
model is even capable of describing the lattice data in the temperature region between Tq and 
1.2 Tc, where the HTL model significantly underestimates the data. 

Finally, figure |5| shows the background field B{T) as a function of temperature. Although 
the setup of the confinement model is quite different from the HTL model, the shape of this 
function remains roughly the same. Note, however, that the B2 term shows a completely 
different temperature behaviour than in the HTL model: it is monotonously decreasing and 
negative from Tc on. The Bi term is vital to reproduce the necessary peak structure. The 
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A 3/4 s/T' 
O 3p/T' 



I ^ ■ ' ' ' ' ' ■ 1 

1 2 3 4 5 

T/Te 

Figure 3: The normalized energy density e = e/T^, entropy density s = 0.75 s/T^ and pressure 
P = 'ip/T^ of our model (solid lines) compared to continuum extrapolated SU(3) lattice data 
(symbols) [Q]. The size of the symbols reflects the lattice uncertainties. The arrow indicates 
the ideal gas limit for massless gluons. 

3 I \ ■ \ ' \ ' \ ■ 




Figure 4: The interaction measure A = (e — 3p) /T^ of the confinement model (solid line) versus 
results (symbols) from a 32^ x 8 lattice. The data symbols represent the continuum interpolated 
values [|]. 
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Figure 5: The background field B(T) and its components Bi and B2, defined in eq.(lS). Also 
shown is B(T) in the HTL model. Symbols display the spacelike plaquette expectation value 
^Ao- taken from the lattice calculation of ref.[Q]. 

maximal value of B{T) is a factor of ~ 1.6 larger than in the HTL model, but it also becomes 
negative for larger temperatures, with its zero Tq slightly shifted from Tq ~ 2 Tc to Tq ~ 2.2 Tc- 
An intriguing observation is that its shape now closely resembles the temperature dependence 
of the spacelike plaquette expectation value Ag-. The space- and timelike plaquettes, Ao- and 
A^, are related to the interaction measure (^) by A = Ao- + A,- and can be expressed in terms 
of the thermal chromomagnetic and chromoelectric condensates {B'^)t and (E^)r as 

^{B')t = -^T'A, + ^{G')t=o and 
vr 11 11 

-(E^)t = ^T'Ar-^{G')T=o. (23) 
vr 11 11 

What we find in fact is B{T) = ^A^{T) (see figure |). This relation between B{T) and (B^)^. 
may be accidental, but it may also hint at a deeper connection between the background field as 
a carrier of non-perturbative effects, and the magnetic condensate. After all, B(T) represents 
the thermal energy of the (non-trivial) Yang-Mills vacuum. 



3 Quasiparticle model with dynamical quarks 

The extension of the mechanism presented in the last section to systems with dynamical quarks 
is not straightforward. Simulations of fermions on the lattice are still plagued by problems. 
No concise continuum extrapolation of the QCD EoS with realistic quark masses exists to 
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date. Nevertheless, it is still possible to construct a model of the EoS with quarks, using some 
reasonable arguments based on the available lattice data. 

3.1 Lattice QCD thermodynamics with quarks 

There have been lattice calculations of the pressure with different numbers of quark flavours 
Nj In the following we focus on results of the Bielefeld group in ref.Q where a p4- 

improved staggered action on a 16'^ x 4 lattice was used. There, the A^-dependence is known 
to be small, in contrast to the standard staggered fermion actions which show substantially 
larger cut-off effects. Lattice calculations were performed for two and three flavour QCD with 
quarks of mass rUq/T = 0.4, and for three flavours with two light quarks [mq/T = 0.4) and an 
additional heavier quark (m^/T = 1.0). From the experience in the pure gauge sector, it has 
been estimated that the continuum EoS lies about 10 - 20% above the data computed on finite 
lattices. 



Figure y displays the lattice pressure, normalized to the Stefan-Boltzmann ideal gas value, 
for the pure gauge system and for systems with 2, 2+1 and 3 quark flavours. A striking 
feature is that, within the errors arising from the cut-off dependence, the QCD EoS shows a 
remarkable flavour independence when plotted against T/Tq- This picture suggests that the 
flavour dependence is well approximated by a term reminiscent of an ideal gas, 

p{T, Nf) oc (^16 + J • 2NcNf'^ ^ p{T/Tc) (24) 

with a function p{T/Tc)- Since Tq changes obviously with the number of degrees of freedom 
present in the thermal system (and therefore with Nf), p is also implicitly A'^^-dependent. 
Scaling against T/Tc, however, the shape of p remains almost the same, indicating that the 
confinement mechanism itself is only weakly flavour-dependent. Therefore, once we understand 
the machinery that is responsible for p{T) in the gluon sector, an extrapolation to systems 
including quarks appears feasible, with the confinement model suitably adapted. Lattice results 
on the order of the phase transition in full QCD support this idea by indicating that the 
transition is first order in the case of three light, degenerate quark flavours and most likely 
second order for two flavours 2C]. Note however that, after applying the previously 

mentioned 10-20% correction to the lattice data, the continuum estimate of the pressure with 
dynamical quarks is much closer to the ideal gas limit than in the pure gauge sector. 



3.2 Thermal masses 

No lattice data on thermal masses with dynamical quarks are available. We thus construct 
effective masses for quarks and gluons by assuming that the Nq- and A^j-dependence of rUg 
and niq are both given by the HTL asymptotic limit. For the thermal gluon mass we employ 
the ansatz: 

^ = lf + §S(T,Nc,N,) (25) 

with the effective coupling 

a(r.A^C,iV;)= ^„^^°_,^^ ([l + i|-^)'. (26) 
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Figure 6: The pressure, normalized to the Stefan-Boltzmann ideal gas value, for the continuum- 
extrapolated pure gauge system and for systems with 2, 2+1 and 3 flavours on a 16^ x 4 lattice, 
obtained with a p4-improved staggered fermion action. The continuum limit is estimated to lie 
about 10 - 20% above the curves shown (flgure from p^). 




go, 6 and /3 are taken to be universal. Setting gQ = 9.4, 6 = 10~^ and (3 = 0.1, the gluon and 
quark masses coincide with the two flavour HTL masses at T ~ 3 Tq, using the parameters of 
ref.Q. The thermal quark mass becomes ^ 



T 



2 



+ \ -3^9{T)\ +-8r^~9{TY (27) 



\\T V 16iVc / 16A^c 



with the zero-temperature bare quark mass rriq^Q. 



In [23 1, a non-perturbative dispersion equation for a thermal quark interacting with the gluon 



condensate has been calculated, and it has been found that the effective quark mass is given 
by ruq ~ 1.15 T in the temperature range between 1.1 Tq and 4 Tc- Eq.(p7[) is within 10% 



in agreement with this result. Nevertheless, very close to Tc the parametrization (|27D may 
be too simple: if the expected chiral phase transition is second order (or weakly first order), 
fermions may decouple in the vicinity of the phase transition because they have no Matsubara 
zero modes, and the transition dynamics would be dominated by the bosonic gluons only. In 
this case gluon masses should become independent of Nj. However, as in the pure gluon sector, 
the results are stable against small variations of the mass parametrizations, and as long as no 
further information is available, eqs.(p7[) and(p^) may be taken as an educated guess. 



We proceed now as follows: First, we assume that the continuum limit of the pressure can 
be obtained from the A^,- = 4 lattice data by applying a 10% correction as mentioned, i.e. 
Pcont — 1-1 Plat- Second, using Occam's razor we employ a umwersa/ confinement function C{T) 
for both quarks and gluons, motivated by eg. (p4|) . The extension of eq.(p!o|) to systems including 



14 





Co 






2 flavours 


1.25 


0.02 


0.28 


2+1 flavours 


1.16 


0.02 


0.29 


3 flavours 


1.03 


0.02 


0.2 


gluon 


1.25 


0.0026 


0.31 



Table 2: Parametrizations of eq.(19) for the confinement function C{T) in the presence of 
dynamical quark flavours. For comparison, the corresponding values of the pure gauge system 
(set B) are also shown. 



quark flavours is straightforward: 



P{T) = ^J dk [C{T)fB{El)] ^ + Y.^j [c{T)fo{ElS\ — - B{T). (28) 



Here = yA;^ + m2(T) as before, and fniE) = {exp{E/T) + 1) ^. The quark energy is 
E^ = Jk'^ + m? for each quark flavour q = i, and mq{T) is given by eq.(p7|) with the bare 



quark masses mq^Q. The background field B{T), the entropy density s and the energy density e 
follow analogously. If the confinement model is applicable, we should expect that the parameters 
of C(r) in eq. ([l^) , as shown in table |l| for the gluonic calculations, are similar in the presence 
of quarks. We start therefore with the gluon values for Co, 5c and (3c and vary them slightly 
until good agreement with the lattice pressure is obtained. To account for the temperature- 
dependent hare masses used specifically in the lattice calculations, the quark masses ruq^ in 
eq.(p7|) are set to mqfi = 0.4 T (light quarks) and rusfi = 1.0 T (heavy quark). Figure |7| shows 
the results for 2 and 2+1 flavours, the corresponding values for the parameters of C(T) can be 
found in table |2|. We observe that indeed, the conflnement factor C(T) does not differ much 
from the pure gluonic case. The factor Bq is set to 1.4 ~ (180 MeV)^. The larger value for 5c 
is explained by noting that, for T < Tq, many light quark-antiquark composites (pions, kaons 
etc.) are present. They contribute sizably to the entropy in the hadronic phase. Accordingly, 
C(Tc) is larger than in the pure gluon case, hence 5c has to increase. 

3.3 Physical quark masses 

Pressure, energy and entropy density for physical quark masses are finally obtained by setting 
niqfi in eq.(p7D to the real- world values mu,d — and ~ 170 MeV. This procedure assumes 
that C(T) is independent of mq^Q which is not clear. In present lattice simulations, the pions 
are too heavy, m^* > 450 MeV, therefore their contribution to s or p is strongly Boltzmann 
suppressed. Since 

g-mi-VTc 1 

^-mi^yyTc ^ 7' 

future computations with lighter, more realistic pion masses are expected to find a substantially 
larger pressure and entropy in the hadronic phase below Tq. With this caveat, we assume for 
now that C(T) does not depend on niq^o. Figure ^ shows a prediction of e, s and p for massless 



15 



5.0 



N,=2+1 




Figure 7: The rescaled lattice pressure pcont — (grey lines) for 2 and 2+1 flavours and the 

pressure obtained from the confinement quasiparticle model with running bare quark masses 
(dashed lines). Values for the parameters of C(T) are shown in table ^. 

two flavour QCD. Reassuringly, the pressure of the confinement model is well within the estimate 
for the continuum EoS of ref.Q for T > 2 Tq- In contrast to the pure gluon EoS, we observe 
that the energy and entropy are close to the ideal gas limit already at T = 3 Tc- However, 
it has to be borne in mind that their normalization is set by Cq which in turn depends on 
the continuum estimate of the N^- = 4 lattice data. More reliable estimates for the continuum 
pressure are needed to confirm this behaviour. It is also worthwhile noting that, going from 
temperature-dependent bare masses (as used in the lattice simulations) to the chiral limit, the 
corresponding change of the pressure in the confinement model is stronger than expected from 
an ideal Fermi gas. It rises by about 7% whereas for an ideal gas with quark mass mq/T = 0.4 
the difference would be only about 3.5% (for Nf = 2). 

In figure ^ we plot the energy density for three light quark flavours on a 16^ X 4 lattice ||1 
and as obtained in the confinement model. Since there is no estimate of the continuum limit 
of these lattice data, the normalization, set by Cq in eq.([l9|), is substantially smaller than in 
previous cases. Apart from that, the data are very well reproduced down to Tc- Finally, figure 



10 shows the calculated pressure, energy and entropy density for a system with two light quark 
fiavours {mq^ = 0) and a heavier strange quark {rngfi — 170 MeV). Here, the approach to the 
Stefan-Boltzmann limit is obviously slower than in the two fiavour case because of the mass 
suppression of the third, heavier fiavour. Nj in eqs.(^), ( p6| ) and ( p7| ) was set to 2.3. The 
results for e and p are also in good agreement for T > 2 Tc with an EoS obtained in the HTL 



quasiparticle model [24|. Closer to Tc, the confinement model predicts a stronger decrease 
of the energy density, though. Again, the background field B{T) resembles the shape of the 
corresponding function in the HTL model. However its zero, Tq, is considerably shifted, from 
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15.0 



N,=2 




1.0 2.0 3.0 

T/T, 

Figure 8: Pressure, energy and entropy density for two light quark flavours in the confinement 
model. The arrow indicates the ideal gas limit. The grey band is an estimate of the continuum 
EoS for massless two flavour QCD, based on an extrapolation of lattice results 

To ~ 1.7 Tc to To ~ 2.7 Tc- 

4 Summary and outlook 

We have presented a novel quasiparticle description of the QCD EoS at finite temperature. 
Our main modification as compared to previous work is the schematic inclusion of confinement 
thas has so far been neglected. First, we have considered thermal SU(3) gauge theory. Guided 
by lattice results for the thermal screening masses, we constructed a thermal mass mg[T) 
for the gluon and calculated the corresponding entropy density as a measure of phase space 
occupation. The difference to the lattice entropy was then attributed to the confinement process 
setting in as Tc is approached from above. Since confinement simply reduces the number of 
thermally active degrees of freedom in a statistical sense, this non-perturbative behaviour is 
incorporated in a model of quasifree, massive quasiparticles by a modification of the particle 
distribution functions with a confinement factor C(T). The energy density and the pressure 
are then uniquely determined. Our model agrees very well with continuum-extrapolated lattice 
data. The interaction measure A(T) is even better reproduced than in the HTL model. As a 
new and interesting aspect, a possible physical connection of the background field B{T) (which 
is the thermal energy of the Yang-Mills vacuum) with the chromomagnetic condensate (B^)^ 
has been proposed. Our results turn out not to be very sensitive to the detailed parametrization 
of the T-dependent quasiparticle mass; the driving new feature is the confinement factor C(T). 

For systems with dynamical quarks, sufficiently precise lattice data are not yet available. With 
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Figure 9: Energy density estimate in the chiral limit for three quark flavours on a 16^ X 4 lattice 
p2|| (data points) and as obtained in the confinement model (solid line). 




Figure 10: Pressure, energy and entropy density for two light quark flavours {ruu^ = 0) and 
a heavier strange quark (m^ ~ 170 MeV) in the confinement model. The arrow indicates the 
ideal gas limit of massless three-flavour QCD. 
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reasonable assumptions, we can nevertheless reproduce existing continuum estimates of lattice 
results for the pressure. Of course, future simulations with higher statistics and smaller quark 
masses are needed to judge if our confinement model yields sensible results in the presence of 
quark flavours. It is encouraging, though, that for temperatures T ^ 2 Tc, the predictions of 
the confinement model with realistic quark masses agree well with HTL model calculations and 
an estimate of the continuum EoS in the chiral limit. 

The function C(T) parametrizes our ignorance about details of the confinement mechanism. It 
would be desirable to connect this macroscopic quantity with microscopic, first-principle QCD 
dynamics, preferably starting in the gluon sector. Interesting questions arise in this context, 
as to the role of instantons and condensates of magnetic monopoles. We point out that a 
comparison of B{T) with lattice data for the spacelike plaquette Aa- in the presence of quarks 
may shed more light on such conjectures. 
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